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Abstract 

An  alternative  and  rigorous  formulation  of  electromagnetic  scattering  — 
the  Delta  Boundary  Operator  (DBO)  technique  —  has  been  reported  in  the 
literature.  Use  of  a  simple  approximation  allows  fast  and  accurate  calcula¬ 
tions  of  scattering  by  rough  planar  surfaces.  Here  an  investigation  of  the 
wider  applicability  of  approximate  DBO  techniques  is  pursued.  Attention  is 
focussed  on  scattering  by  simple  smooth  bodies  in  the  so-called  resonance  re¬ 
gion,  where  the  approximations  of  high  frequency  asymptotic  techniques  are 
often  inadequate.  It  is  shown  that  application  of  approximate  DBO  meth¬ 
ods  leads  to  predictions  of  non- vanishing  surface  currents  in  unlit  regions  and 
smooth  transition  currents  at  shadow'  boundaries. 


Copyright 

C 

Controller 
HMSO  London 


1 


INTENTIONALLY  BLANK 


Contents 

1  Introduction 


1 


2  Principles  of  the  Rigorous  DBO  Technique  2 

2.1  Scattering  by  an  infinite  perfectly  conducting  plane .  4 

2.2  Scattering  by  a  perfectly  conducting  circular  cylinder .  5 

3  Principles  of  the  Approximate  DBO  (ADBO)  Technique  7 

4  Application  of  ADBO  Theory  to  Scattering  by  a  Circular  Cylinder  10 

5  Analysis  of  the  r/’-Distribution  of  the  Circular  Cylinder  14 

6  Application  of  IDBO  Theory  to  Scattering  by  a  Circular  Cylinder  18 


7  Prescription  of  the  IBBO  Technique  for  a  General  Two-Dimensional 


Smooth  Scatterer  21 

8  Conclusions  and  Suggestions  for  Further  Work  23 

9  References  25 

A  Analysis  of  the  Integral  27 

A.l  Stationary  Phase  Point  Contributions .  28 

A. 2  End-Point  Contributions .  30 

B  Asymptotic  {ka  1)  Analysis  of  t/’‘^(s,0)  33 

C  Analysis  of  the  Integral  /:a)  37 

C.l  Saddle  Point  Contributions .  37 

C.2  End-Point  Contributions .  38 

D  Contribution  of  to  the  Cylinder  Current  39 


Accession  For 


INTENTIONALLY  BLANK 


1  Introduction 


Electromagnetic  scattering  by  bodies  up  to  a  few  wavelengths  in  size  can  be  calcu¬ 
lated  using  wire  grid  or  surface  patch  modelling  in  conjunction  with  the  method  of 
moments.  High  frequency  techniques  such  as  the  Geometrical  Theory  of  Diffraction 
(GTD)  or  Physical  Optics  (PO)  supplemented  by  a  form  of  the  Physical  Theory  of 
Diffraction  (PTD)  can  often  be  successfully  applied  to  larger  objects.  However  with 
complex  scatterers  there  is  often  an  intermediate  region  where  the  computational 
requirements  of  moment  method  calculations  are  too  great  and  where  the  approxi¬ 
mations  of  the  high  frequency  asymptotic  techniques  are  inadequate.  Scattering  by 
bodies  of  this  size  is  a  topic  currently  attracting  a  great  deal  of  interest. 

An  alternative  and  rigorous  formulation  of  electromagnetic  scattering  -  called 
the  Delta  Boundary  Operator  (DBO)  technique  -  has  been  reported  by  Maystre  [1] 
and  by  Saillard,  Roger  and  Maystre  [2].  Central  to  this  scheme  is  the  solution  of  an 
abstract  boundary  value  problem  from  which  the  solution  of  a  physical  scattering 
problem  can  be  derived.  It  has  been  shown  [3,4]  that  a  simple  approximation  to 
the  solution  of  the  interm.ediate  problem  allows  fast  and  accurate  calculations  of 
scattering  by  rough  planar  surfaces  where  the  scale  of  the  roughness  is  of  the  order 
of  the  wavelength  of  the  incident  radiation.  In  particular  the  solutions  are  more 
accurate  than  those  obtained  by  using  Beckmann’s  method  [5,6],  a  high  frequency 
technique  which  uses  the  Kirchoff  approximation.  In  this  paper  an  investigation 
of  the  wider  applicability  of  approximate  DBO  techniques  is  pursued.  Attention  is 
focussed  on  scattering  by  simple  smooth  bodies.  Emphasis  is  given  to  scattering  by 
objects  in  the  resonance  region  and  comparison  is  made  with  the  PO  technique. 

In  Section  2  the  principles  of  the  rigorous  DBO  formulation  are  outlined  and 
illustrated  with  simple  examples.  A  simple  approximate  DBO  (ADBO)  theory  is 
developed  in  Section  3  and  this  is  applied  to  plane  wave  scattering  by  a  conducting 
circular  cylinder  in  Section  4.  Some  refinement  is  seen  to  be  necessary  for  accurate 
current  and  radar  cross  section  (RCS)  predictions  and  an  improved  DBO  (IDBO) 
procedure  is  proposed  in  Section  5.  This  is  tested  by  apphcation  to  the  circular 
cylinder  problem  in  Section  6  and  is  prescribed  for  a  general  smooth  two-dimensional 
{2D)  scatterer  in  Section  7.  Section  8  includes  suggestions  for  further  validation  and 
development  of  the  IDBO  scheme. 
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2  Principles  of  the  Rigorous  DBO  Technique 

This  section  draws  heavily  on  the  work  of  Maystre  [1].  It  is  included  to  establish 
notation  and  to  fatniharise  the  reader  with  the  nature  of  the  method. 

Consider  the  scattering  of  a  plane  electromagnetic  wave  incident  normally  upon 
an  infinitely  long  cyhndrical  body  situated  in  free  space,  as  shown  in  Fig.l.  The 
2-axis  is  chosen  to  be  parallel  to  the  axis  of  the  cylinder.  A  general  point  on  the 
surface  contour  C  will  have  arc  length  s  from  a  reference  point  P. 

The  body  will  be  assumed  to  be  perfectly  conducting.  The  incident  radiation 
is  assumed  time-harmonic.  All  fields  will  then  exhibit  the  same  time  dependence, 
chosen  to  be  which  will  be  suppressed  throughout.  We  confine  the  analysis  to 

the  case  where  the  incident  electric  field  is  2-polarised.  (This  case  is  also  commonly 
referred  to  as  ^'-polarisation  or  TM-polarisation.)  Then  the  total  electric  field  ev¬ 
erywhere  has  a  non-zero  2-component  only.  At  a  point  external  to  the  scatterer  with 
position  vector  r  this  is  given  by 

£(r)  =  E'(r)  +  E‘(r),  (1) 

where  E'(r)  and  ^^(r)  denote  the  incident  and  scattered  fields  respectively.  The 
scattered  field  satisfies  the  Helmholtz  equation 

(V"  +  l:")£'(r)  =  0,  (2) 

where  and  /xq  and  Cq  are  respectively  the  permeability  and  permittivity 

of  free  space.  In  addition  the  scattered  field  is  subject  to  the  boundary  condition 

Enc=-E'\c  (3) 

on  the  surface  C  and  to  the  far-zone  Sommerfeld  radiation  condition. 

The  above  completely  specifies  a  physical  electromagnetic  scattering  problem 
and  in  principle  a  unique  solution  for  the  electric  field  may  be  obtained.  Often  the 
surface  current  density  J{s)  is  determined  first.  The  scattered  field  outside  the  body 
is  then  calculated  by  means  of  the  Helmholtz  representation  : 

E^(r)  =  -  [  G(r,r');(r')ds',  (4) 

Jc 

where  r'  =  t'{s')  is  the  position  vector  of  a  general  point  on  the  contour  C,  j  = 
—iujfioJ  represents  a  normalised  current  density  and  the  two-dimensional  Green’s 
function  is  given  by 

G(r,!;')=  (5) 

Now  consider  the  following  abstract  problem  —  called  the  delta  boundary  opera¬ 

tor  (DBO)  problem  by  Maystre.  A  fictitious  source  outside  the  scatterer  is  assumed 
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to  be  such  that  it  gives  rise  to  an  incident  field  U*  of  unit  amplitude  at  a  point  s  =  Sq 
on  the  surface  of  the  body  and  which  vanishes  elsewhere  on  C;  i.e. 

U'  Ic  =  -  So).  (6) 

(See  Fig.2.)  The  presence  of  the  scattering  object  causes  a  scattered  field  U*  w’hich 
satisfies  the  following  conditions  : 

(V^  +  k‘^)W{r)  =  0  outside  C  (7) 

U>\c  ^-8{s-so).  (8) 

In  addition  the  scattered  far-held  must  obey  the  appropriate  radiation  condition. 
The  DBO  problem,  as  stated  here,  is  simply  a  mathematical  boundary  value  prob¬ 
lem.  A  unique  solution  for  the  scattered  field  is  possible  in  principle. 

Maystre  defines  a  so-called  ^’-distribution  for  the  DBO  problem  which  is  analo¬ 
gous  to  the  current  j  of  the  physical  problem; 

dU‘ 

V'(s,so)  = (9) 

The  solution  of  the  physical  scattering  problem  may  be  obtained  from  that  of  the 
DBO  problem  by  means  of  the  prinrq.h  *■  superposition.  The  required  superposition 
is  found  by  inspection  of  the  bounda:  conditions  (3)  and  (8)  and  by  use  of  the 
relation 

E'{s)=^  J^6{s-so)E'{so)dso.  (10) 

For  example,  the  scattered  field  on  the  surface  of  the  scatterer  is  given  by 

E’(s)  =  [  U’{s,so)E'{so)dso  =  -E'{s)  (11) 

Jc 

as  required.  Similarly  the  current  j{s)  is  given  by 

8F'  r 

=  J^xl^{s,so)E'{so)dso.  (12) 

Thus  the  physical  scattering  problem  may  be  solved  by  first  solving  the  associated 
DBO  problem  to  obtain  the  ^’-distribution. 

It  is  important  to  appreciate  that  the  rigorous  DBO  method  of  solution  of  a  scat¬ 
tering  problem  offers  nothing  more  than  an  alternative  analysis.  It  is  not  anticipated 
that  application  of  the  method  will  lead  to  new  exact  analytic  results  for  previously 
unsolved  problems.  As  will  become  clear  in  Section  3  the  value  of  this  approach 
is  that  it  suggests  approximate  techniques  which  are  expected  to  be  both  fast  and 
accurate  for  many  scattering  problems.  With  this  in  mind  we  now  demonstrate  the 
DBO  method  by  applying  it  to  two  simple  problems  with  known  solutions. 


2.1  Scattering  by  an  infinite  perfectly  conducting  plane 

The  plane  is  taken  to  be  coincident  with  the  y  =  0  plane.  We  consider  radiation  in 
the  region  y  >  0.  The  DBO  problem  is  specified  as  follows  : 

{W^  +  P)U%x,y)  =  0  fory>0  (13) 

U‘ix,0)  =  -6{x-xo)  (14) 

U’{x,y)  satisfies  a  radiation  condition  as  y  — »  +oc-.  (15) 

To  obtain  a  solution  we  express  U*  as  a  Fourier  integral  : 


U%x,y)=  r  Viy.hy^^dh. 

J  -'OC 


Substitution  into  the  Helmholtz  equation  gives 

(PV 

+  {P  -  y)V  =  0  for  y  >  0. 

ay2 

The  solution  for  V’  which  ensures  compliance  with  the  radiation  condition  is 

T(y,A)  =  V(0,/l)e-^‘'••*)^ 


By  using  the  relation 


for|/.|>A- 

j  i0T\k\<k. 

s{x-xo)  =  ~  r 


imposition  of  the  boundary  condition  (14)  gives 


Thus  we  deduce 


Q.7r  J —(x> 


In  particular  the  i/’-di^tribution  for  this  problem  (denoted  is  given  by 

r{x,Xo)  =  V(x,0)  =  —  /  (23) 

uy  Z'K  y  — oo 

Manipulations  with  the  Fourier  transforms  of  Hankel  functions  allow  closed  form 
evaluation  of  this  integral  [1,7].  We  find 

=  (2.) 

2  (x  -  xo) 
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Observe  that  0p(x,io)  depends  only  on  x  —  xq,  a  consequence  of  the  translation 
invariance  of  the  scatterer.  For  a  general  cylindrical  body  Maystre  [1]  has  shown 
that  the  ^-distribution  is  a  symmetric  function;  i.e. 

xp{s,So)  =  ip{so,s).  (25) 


This  is  a  manifestation  of  the  reciprocity  principle. 

Now  consider  the  problem  of  scattering  of  a  plane  wave  incident  in  the  x  —  y 
plane  at  an  angle  q  to  the  surface  normal.  We  have 

£’(x,j/j  =  (26) 


Then  the  current  is  found  from  (‘2)  to  be 


j(^) 


-il-cosQ  - 

-2ikcosa 


/  .  C  uXq 

-  ^o) 


(27) 


in  accord  with  the  well-known  solution.  Here  use  has  been  made  of  the  Fourier 
transform  of  the  Hankel  function  ; 


(28) 


2.2  Scattering  by  a  perfectly  conducting  circular  cylinder 

The  circular  symmetry  of  this  scatterer  implies  a  t/>-distribution  V’'^(s,so)  dependent 
only  on  s  —  Sq.  Without  loss  of  generality  we  may  thus  set  sq  to  zero  and  determine 
0). 

For  a  cylinder  of  radius  a,  s  =  a9  and  the  DBO  problem  may  therefore  be  stated 
as  follows  : 

(V^  +  P)U%r,  e)  =  0  for  r  >  a  (29) 

U^{a,e)  =  -6{e)/a  (30) 

U’{r,B)  satisfies  a  radiation  condition  as  r  — +  oc.  (31) 

A  straightforward  calculation  [1]  produces  the  result 


r{s,0) 


n^oc  Hi!\ha) 


(32) 


For  a  plane  wave  incident  from  the  -f  x-direction  (see  Fig. 3) 


(33) 
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and  (12)  gives  the  surface  current 


j{ae)  =  ^{a,  e)  +  a  aeo)E'{a,  0o)d0o- 

or  Jo 

By  expressing  E'  as  a  Fourier  series  : 

E'{r,e)=  £  (-t)V„(A:r)e-"^ 

ni--oo 

it  can  be  readily  verified  that 


.  2i  ^ 

;(««)  = -j;  E 

n=-oc 


Thus  the  classical  solution  is  reproduced.  (See,  for  example.  Ref.8.) 


3  Principles  of  the  Approximate  DBO  (ADBO) 
Technique 

Maystre  [1]  has  performed  analjlical  and  numerical  studies  of  the  functions  rh‘^{s,0) 
and  t/’P(s,0)  =  {—ik/2s)H[^\ks).  His  main  observations  are  listed  below  : 

1.  and  0^(5, 0)  exhibit  identical  singularities  as  s  — >  0.  (Note  however 

that  is  a  27ra-periodic  function  and  so  is  singular  as  s  — >  2m7ra,  where 

m  is  an  integer.  We  limit  comparison  to  the  region  A-jsl  <  irka.) 

2.  V’‘^('S,0)  — ♦  xl>^{s,0)  as  fca  — >  oo  . 

3.  Provided  ka  ~  2,  v'’‘^(s,0)  «  vP(s,0)  for  /:|sj  ~  nka. 

4.  Values  of  both  the  real  and  imaginary  parts  of  become  negligible  for 

^'|s|  ~  27r. 

5.  Values  of  both  the  real  and  imaginary  parts  of  become  negligible  for 

k\s\  ~  27r  provided  ka  ~  2. 

Recalling  the  relationship  (12)  between  the  ^’-distribution  and  the  surface  current, 
observations  4  and  5  suggest  the  following  fundamental  conclusion  [1,3]  for  TM 
incidence  on  a  general  2D  scatterer:  The  surface  current  at  any  point  on  the  body 
is  influenced  by  the  shape  of  the  scatterer  within  a  distance  along  the  surface  of 
approximately  a  wavelength.  This  phenomenon  of  short-range  coupling  is  excluded 
from  the  more  popular  Physical  Optics  (PO)  theory  which  gives  rise  to  surface 
currents  which  depend  only  on  the  local  slope  of  the  surface. 

The  observed  similarities  between  the  distributions  v>^(s,0)  and  ^’‘^(s,0)  (for 
ka  ^  2)  suggest  an  approximate  analysis  of  scattering  by  a  general  2D  body  which 
accounts  for  the  effects  of  short-range  coupling.  For  a  scatterer  of  perimeter  I  it  seems 
reasonable  to  expect  (and  assume)  that  the  true  i/’-distribution  i/’(s,so)  (which  will 
be  an  /-periodic  function  in  both  s  and  so)  can  be  accurately  represented  by  the  func¬ 
tion  t/’’’(s,  So)  when  |s  —  sqI  <  1/2.  This  we  shall  call  the  approximate  DBO  (ADBO) 
approach.  Underlying  this  assumption  are  two  quite  separate  approximations: 

1.  It  is  assumed  that  xl){s,So)  depends  on  the  variables  s  and  So  solely  through 
their  difference  s  —  So;  i.e.  we  assume  r/’(s,5o)  ~  t/’(s  —  so,0).  This  is  certainly 
justifiable  for  bodies  with  surfaces  which  are  small  perturbations  on  a  plane  or 
circular  cylinder  but  the  validity  of  its  wider  applicability  is  uncertain. 

2.  It  is  assumed  that  ih(s  —  soiO)  «  V’^(s  —  so<0)  =  T/)p(s,so)  for  js  --  sqI  ^  //2. 
In  view  of  observations  1  and  3  on  the  ^’-distribution  of  the  circular  cylinder, 
it  seems  reasonable  to  anticipate  that  this  approximation  will  introduce  little 
error  provided  the  scatterer  is  smooth  (i.e.  edges  are  absent)  and  is  not  too 
small;  i.e.  /  ~  2A  or  /:/  ~  47r.  (A  =  2n/k  denotes  the  wavelength  of  the  exciting 
radiation.) 


The  assumption  of  a  symmetric  ^-distribution  ensures  that  fields  obtained  using 
this  approximate  technique  obey  the  reciprocity  principle.  (See  Section  2.)  It  is 
worth  noting,  though,  that  approximate  solutions  satisfying  reciprocity  may  not  be 
any  more  accurate  than  non-reciprocal  ones  (9]. 

Since  xP{3,sq)  is  /-periodic  in  sq  (and  s),  (12)  gives 


.  dE' 

j{s)  =  -^{s)  +  /  xlf{s,SQ)E'{so)dso. 
on  Jm-ii2 


rt+t/2 

-1/2 


Now,  invoking  approximation  1  above,  we  write 

dE'  r‘+‘n 


OE' 

j{s)  ^  rl’is  -  so,0)E'{so)dso 

an  Js-i/2 


8  F*  r'/2 

—  (s)-h  /  rlis\0)E'is'-\-s)ds'. 

on  J-U2 


-1/2 

r'/2 


dn  '  J-i/2 

Implementation  of  approximation  2  then  gives 


(37) 


(38) 


d  F'  FF 

j{s)  ss  jADBo{s)  =  -^(s)  +  /  rl’^{s',0)E'{s'  -1-  s)ds'. 

dn  J-i/2 


(39) 


Recalling  that  t/)P(s,0)  is  numerically  negligible  for  |s|  ~  A  (observation  4)  we  see 
uhat  little  loss  of  accuracy  will  ensue  by  changing  the  integration  limits  to  ±oo;  i.e. 

d  F'  r°^ 

jADBois)  «  +  J_^rl’ns\0)E'{s'^s)ds'.  (40) 


This  further  approximation  is  expected  to  be  valid  provided  /  ~  2A. 

The  replacement  of  the  true  V’-distribution  by  its  counterpart  for  the  plane  may 
be  viewed  as  an  infinite  plane  approximation.  PO  uses  an  infinite  tangent  plane 
approximation  to  determine  the  current  locally  and  the  shape  of  the  scatterer  is 
accounted  for  when  integrating  the  current  to  find  the  RCS.  However,  in  the  ADBO 
technique  the  infinite  plane  approximation  is  used  one  stage  earlier,  allowing  the 
current  to  be  influenced  by  the  geometry  of  the  scatterer  in  the  nearby  region. 

It  seems  reasonable  to  expect  (40)  to  be  an  accurate  representation  of  reduced 
surface  current  on  many  2D  scatterers.  Furthermore  the  current  is  expressed  as 
a  convolution-type  integral  which  can  be  performed  using  Fourier  transform  tech¬ 
niques.  Rapid  calculation  of  currents  can  therefore  be  achieved  by  means  of  fast 
Fourier  transform  (FFT)  algorithms.  Since  the  Fourier  transform  of  is  calcula¬ 
ble  analytically  (see  (28))  this  procedure  also  alleviates  problems  which  w'ould  be 
encountered  in  the  direct  numerical  evaluation  of  the  integral  due  to  the  singular 
behaviour  of 

In  Refs. 3  and  4  this  approximate  technique  has  been  applied  to  scattering  by 
planar  rough  surfaces.  In  the  resonance  region,  where  the  roughness  dimensions  and 
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the  wavelength  of  the  incident  radiation  are  of  the  same  order  of  magnitude,  the 
results  were  consistently  more  accurate  than  those  obtained  by  using  the  Beckmann 
method  [5,6],  which  strictly  is  applicable  only  when  the  wavelength  is  small  relative 
to  the  scale  of  roughness.  In  Section  4  the  application  of  the  ADBO  method  to  a 
simple  smooth  scatterer  —  the  circular  cylinder  —  is  investigated  in  detail. 
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4  Application  of  ADBO  Theory  to  Scattering  by 
a  Circular  Cylinder 

The  scattering  of  a  TM-polarised  plane  wave  by  a  perfectly  conducting  circular 
cylinder  has  been  examined  rigorously  in  Section  2.  Here  the  problem  is  revisited 
using  the  ADBO  theory.  From  the  arguments  of  the  previous  section  we  expect 
reasonable  accuracy  for  ka  ~  2.  From  (39)  we  find 


jADBoi^)  =  —ik  cos  6e 


^ika  co£  $ 


ik  ^  ^  ,  . 

-  y 


As  a  first  attempt  at  an  analysis  of  the  current  we  invoke  the  additional  approxi¬ 
mation  7,(0,  ia)  »2  Iao{6,ka)  in  accordance  with  (40).  For  an  approximate  analytic 
evaluation  of  the  integral  loo  consider  the  following  argument.  From  observation  4 
of  Section  3  we  see  that  the  distribution  H\^\ka6')  j 6'  is  numerically  significant  only 
when  \6'\  ~  2-k fka.  In  view  of  this  we  are  able  to  replace  by  its  small 

|0'|  approximation  (strictly  valid  for  1^'|  1)  when  ka  27r.  Thus 


/oo(0,  ka) 


^-ikacose  ^\kae' sine  jni 


=  2|cos01e 


cos  B 


where  use  has  been  made  of  the  standard  integral  (28).  Therefore,  when  ka  ^  2it, 


jADBoi^) 


-2ziccos0e-*“'°*'’  for  0  <  |01  <  7r/2 
0  for  7r/2  <  1^1  <  TT 


These  expressions  for  the  reduced  current  on  a  conducting  cylinder  are  precisely 
those  predicted  by  PO  theory.  However,  a  more  rigorous  analysis  of  /„  is  possible 
and  is  given  in  Appendix  A.  There  it  is  found  that 


h{0,ka) 


ka^\,cos$^0 


21  cos  0|e 


— lito  COS  9 


2ka  cos^  6 , 


where 


2  \3/2  jA^inka  ^ikacosB 

ka)  7r^cos^0 


1/2  /  0  \ 


giA:a(3r/2-|91)  g.fca(!e|-r/2)  \ 

"(37r/2  -  101)3/2  (101  -  7r/2)3/V  ■ 
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U{x)  denotes  the  unit  step  function  and  Ai{x)  represents  the  Airy  function.  There¬ 
fore  the  ADBO  theory  gives  the  following  expressions  for  the  current  on  a  circular 
cyhnder: 


-2.A:cosee-‘“'“*  (l  + 

-'inieM) 


for  0  <  |0|  ~  7r/2 

for  Tl  12  ~  1^1  <  TT 


(48) 


In  accordance  with  (44)  PO  currents  are  recovered  in  the  high  frequency  {ka  — ►  oo) 
limit.  These  results  are  valid  provided  ka  >>  1  and  cos  9  ^  0.  Therefore  we  have 
been  unable  to  derive  a  simple  analytic  form  for  the  current  in  the  transition  region 
|0|  Rs  7r/2.  Outside  the  shadow  boundary  region  we  have  obtained  expressions  for 
correction  terms  to  the  PO  currents  which  account  for  the  effects  of  the  nearby 
scatterer  shape.  It  is  of  interest  to  compare  (48)  with  the  corresponding  asymptotic 
form  of  the  exact  current  jEx{6)  [10]: 


jEx{6) 


ka'yi  .cot  BijtO 


-2i«:cos«e-‘““‘*  (l  +  5j^) 

forO<|»|iV2 


I  (^)''' j:”-.  ^  i«i  < » 

(49) 

The  infinite  series  represent  creeping  wave  contributions.  We  see  that  the  curvature 
correction  term  in  the  illuminated  region  predicted  by  the  ADBO  method  is  of 
the  correct  form  and  so  breaks  down  in  similar  circumstances  to  jsxi^)-  However 
its  magnitude  is  half  the  true  value.  Note  also  that  l^^{0,ka)  clearly  represents  a 
creeping  wave  contribution  to  the  ADBO  current. 


for  7r/2 


To  check  the  validity  of  the  analytic  approximation  (48)  away  from  the  vicinity 
of  the  shadow  boundary  a  numerical  evaluation  of  /,  is  necessary.  This  can  be 
accomplished  only  after  careful  treatment  of  the  singularity  at  6'  =  0.  The  details 
are  discussed  in  Appendix  A.  I.^{6,ka)  has  been  evaluated  using  (83).  Alternatively, 
and  as  discussed  at  the  beginning  of  this  section,  calculation  of  Ir{9,  ka)  can  proceed 
approximately  through  evaluation  of  loo{0,ka).  (It  is  shown  in  Appendix  A.l  that 
the  difference  /oo(^,  ka)  —  /,(0,  ka)  is  of  order  0(l/fca)^/^  and  therefore  small  relative 
to  either  loo  or  when  fca  >  1.)  The  integral  loo  has  been  evaluated  numerically 
in  three  different  ways: 
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1.  Through  computation  of  the  expression  equivalent  to  (83)  for  lA.{0,ka): 

lA{e,ka)=  -  + 

-  iH^^\kaA){sm{A  +  +  sin(yl  - 

+  i  r  cos(0'  +  0)(1  -  ika  cos(0'  +  e))de'. 

(50) 


Values  of  I^a  have  been  obtained  by  examining  for  convergence  for  large  A. 

2.  By  using  Fourier  transform  techniques  (since  Ioo{0,  ka)  is  a  convolution  type 
integral).  The  equivalent  results  are  obtained  directly  by  means  of  Fourier 
series.  Upon  using  the  identity  (35)  we  find 


loo{0,ka) 


n=-oo  ^ 


de' 


(51) 


where  the  standard  integral  (28)  has  been  used.  Summation  of  this  series  has 
been  performed  by  rewriting  it  as 

2i  °° 

loo{0i  ka)  =  —  ^n(— ka)Jn{ka)  cos  n0,  (52) 

(co  =  li  =  ^2  =  '  •  •  =  2)  and  retaining  sufficient  terms  in  the  sum  to  ensure 
convergence. 

3.  By  using  an  FFT  algorithm. 

The  results  obtained  by  these  three  methods  differ  negligibly. 

Plots  of  the  variation  of  the  magnitude  and  phase  of  the  surface  current  with  angle 
6,  as  given  by  (41),  are  shown  in  Fig. 4  for  ka  =  5  and  =  10.  Also  shown  are  the 
results  with  /,  replaced  by  loo  and  the  analytic  form  (48)  of  the  currents.  As  expected 
invoking  the  extra  approximation  /,  %  loo  introduces  little  error  (particularly  for 
large  ka).  It  is  clear  also  that  (48)  is  indeed  a  valid  representation  of  the  ADBO 
cylinder  current  in  the  deep  illuminated  and  deep  shadow  regions. 

In  Fig. 5  the  magnitude  and  phase  of  the  cylinder  current  obtained  using  the 
ADBO  technique  (by  numerical  evaluation  of  7^)  are  compared  with  the  correspond¬ 
ing  exact  [10]  and  PO  solutions.  In  general  the  ADBO  currents  compare  favourably 
with  the  PO  results.  Unlike  PO,  application  of  the  ADBO  method  leads  to  pre¬ 
dictions  of  non-vanishing  surface  currents  in  unlit  regions  and  smooth  transition 
currents  at  shadow  boundaries.  However  the  ADBO  currents  display  an  erroneous 
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oscillatory  variation  with  angle  for  6  ~  7r/2  and  their  magnitudes  are  over-estimated 
in  the  deep  shadow  region. 

We  now  turn  to  a  determination  of  the  scattered  far-held  (and  hence  the  radar 
cross  section  \^RCS))  using  the  ADBO  technique.  This  is  achieved  by  integration  of 
the  current  around  the  contour  of  the  scatterer  as  indicated  in  (4).  For  the  case  of 
the  circular  cyhnder  the  scattered  far-held  is  found  to  be 

y/8ir  kr  Jo 

Since  we  have  been  unable  to  derive  a  reliable  analytic  expression  for  the  ADBO 
current  at  all  points  on  the  surface  we  revert  to  the  form  (41).  Inserting  into  (53) 
and  interchanging  the  order  of  integration  in  the  second  term  gives 

E\r,B)  -Z!^e‘<*^^-^/‘‘H-2rcos(^/2)Ji(2/l:a  cos{e/2))  +  L40,ka)},  (54) 

ySirkr 


where 


A  H\^\kae') 


6' 


Jo  2ka  cos 


de'. 


The  bistatic  RCS  cr{0)  is  dehned  as 


(y(0)  =  lim  27rr 

r— *00 


E’{r,e) 

E'{r,e) 


(55) 


(56) 


ADBO  predictions  of  bistatic  RCS  are  shown  in  Fig. 6  for  cylinders  of  size  ka  =  5 
and  fca  =  10.  Also  shown  are  the  corresponding  exact  [8]  and  PO  results.  As  with 
care  has  been  taken  with  the  singularity  of  the  integrand  for  the  numerical  evalu¬ 
ation  of  Tt-  There  is  evidence  that  the  over-estimated  currents  in  the  deep  shadow 
region  contribute  signihcantly  to  the  RCS,  with  their  erroneous  oscillatory  angu¬ 
lar  behaviour  persisting  after  integration.  It  is  debatable  whether  application  of  the 
ADBO  technique  to  the  circular  cyhnder  scattering  problem  has  led  to  improvements 
over  the  PO  RCS  predictions. 

Ill  view  of  this  disappointing  conclusion  it  is  natural  to  ask  how  the  ADBO 
scheme  might  be  modified  .o  improve  its  accuracy.  The  two  basic  assumptions 
inherent  in  the  technique  were  given  in  Section  3.  It  has  already  been  remarked  that 
the  i/:- distribution  for  the  cylinder  xl!‘^{s,so)  depends  only  on  s  —  so  ;  i.e.  xl!'^{s^so)  = 
—  SoiO)-  Thus  only  the  second  approximation  was  used  in  the  cylinder  analysis; 
i.e.  we  assumed  ^’‘^(•5,0)  ^  V'^(-s,0)  for  |s|  <  tto.  Although  observation  3  of  Section 
3  supports  this  assumption  it  is  clear  that  it  is  inadequate  for  the  Ccises  we  have 
considered  {ka  =  5  and  ka  =  10).  It  seems  likely  that  a  better  approximation  to 
t/’‘^(s,0)  will  lead  to  improved  current  and  RCS  predictions.  Therefore  in  Section  5 
we  examine  in  detail  the  i/j-distribution  of  the  circular  cylinder. 
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5  Analysis  of  the  ^’-Distribution  of  the  Circular 
Cylinder 


In  this  section  we  seek  an  improved  approximation  for  0).  Since  ^^*^(5,0)  — ♦ 
V’’’(s,0)  as  Aro  — ♦  oo  (observation  2  of  Section  3)  we  look  for  correction  terms  by 
examining  the  asymptotic  behaviour  of  the  exact  expression  (32).  By  using  a  Watson 
transformation  (see  Appendix  B)  we  find 


0)  "  « '  V{s,  0)  -  0), 

4a 


ik 


(57) 


where 

^  n=l  m  =  l 

with 


(58) 

(59) 


0)  represents  creeping  wave  contributions.  For  each  value  of  n  (corresponding 
to  a  root  i/n  of  Hl^J{ka)  =  0)  there  is  an  infinite  number  of  such  terms,  each  integral 
change  in  m  corresponding  to  a  complete  circulation  of  the  cylinder.  However, 
careful  inspection  shows  that  not  all  of  the  expected  creeping  wave  terms  are  present. 
Recalling  that  |0|  <  tt  we  see  in  particular  that  all  waves  emanating  from  the  point 
of  excitation  on  the  surface  and  traversing  less  than  half  the  cylinder  circumference 
are  excluded.  Indeed  ^’“^^^’(■5,0)  offers  a  complete  representation  of  creeping  waves 
only  when  6  =  0.  The  imaginary  parts  of  the  determine  the  rates  at  which 
the  creeping  waves  are  attenuated  as  they  move  around  the  cylinder.  From  (135), 
/m(i/„)  oc  (Aa)’/^  and  so  for  large  ka  there  is  considerable  attenuation.  Bearing 
in  mind  the  remarks  made  above  it  follows  that  the  most  significant  creeping  wave 
contributions  are  excluded  from  Indeed,  even  for  moderate  values  of  ka 

it  is  readily  verified  using  (135)  and  tables  of  the  constants  q„  (see,  for  example. 
Refill) that 

|e"'"(2mi^+®)|,  |g.V„{2mK-e)j  <  2Q-4 


for  all  allowed  values  of  m  and  n.  It  is  immediately  apparent  that  the  contribution 
of  V’‘^^^^(5,0)  to  t/>‘^(s,0)  is  relatively  small  and  indeed  numerically  negligible.  Thus 
in  our  numerical  investigation  of  (57)  we  will  ignore  the  contribution  of  0). 

We  must  however  anticipate  that  the  issue  of  the  ‘missing’  creeping  waves  is  an 
important  one  since  their  numerical  contributions  may  be  significant.  This  matter 
wiU  be  discussed  again  later  in  this  section. 

The  validity  of  (57)  and  the  significance  of  the  curvature  correction  term  can  be 
determined  by  numerical  comparison  with  the  exact  form  (32)  and  the  plane  approx¬ 
imation  t/>^(s,0).  Maystre  [1]  has  rewritten  (32)  in  a  form  suitable  for  computation: 
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= 


2{ka)^sm^ie/2) 


+  In 2 cos  0  sin  |0| 


+  cos^ln  I  sin(0/2)j + 


1  Hi^\ka)  2  Hi^\ka) 
kaH^^\ka)  kaHl^){ka) 


cos  6 


(61) 


The  infinite  series  converges  rapidly.  As  before  6  =  s/a  and  |0|  <  tt.  Plots  of  the  real 
and  imaginary  parts  of  V’‘^(s,0)  are  given  in  Fig. 7  for  the  cases  ka  =  5  and  ka  =  10. 
Also  shown  is  the  plane  approximation  and  the  approximation  to  i/i‘'(s,0) 

given  by  (57).  Two  remarks  are  appropriate: 


1.  The  discrepancies  between  ^^(5,0)  and  the  exact  form  (61)  of  V)‘^(s,0)  are  suffi¬ 
ciently  large  to  support  our  conclusion  that  the  ADBO  procedure  is  inadequate 
for  accurate  current  and  RCS  predictions. 


2.  The  inclusion  of  the  curvature  correction  leads  to  a  significant  improvement  to 
our  approximation  to  t/i‘^(s,0)  only  for  k\s\  ~  1. 

It  is  apparent  from  Fig. 7  that  the  representation  (57)  of  V>^(s,0)  fails  for  fc|s|  ~  1. 
This  is  because  the  magnitude  of  the  curvature  correction  does  not  decrease  with 
increasing  k\s\.  It  is  this  term  which  dominates  for  k\s\  ~  1  giving  rise  to  the 
observed  oscillatory  behaviour.  Clearly  for  a  more  accurate  representation  of  0) 
this  correction  term  must  be  strongly  damped  as  t|s|  is  increased.  The  origin  of 
the  absent  attenuation  is  alluded  to  in  Appendix  B.  There  it  is  noted  that  the 
approximation  (132)  to  /(i/,  ka)  =  Hl^^'{ka)/ Hl^^{ka)  (used  to  determine  (57))  does 
not  have  the  correct  pole  structure.  Nevertheless  (132)  is  a  good  approximation  to 
f{i/,ka)  provided  |i/|  76  ka.  In  particular  it  is  valid  when  |i^|  ka.  That  this  is  so 
is  borne  out  by  the  accuracy  of  the  representation  (57)  of  t/'‘^(s,0)  for  k\s\  ~  1,  since 
the  variables  ^  and  ks  form  a  Fourier  transform  pair. 

To  obtain  a  more  accurate  model  of  ^>‘^(■5,0)  for  k\s\  ~  1  we  must  invoke  a 
representation  for  f{v,ka)  which  has  the  correct  behaviour  when  |i/l  ~  ka.  In 
particular  we  require  that  f{v,  ka)  has  poles  at  1/  =  ±z/„  (since,  by  (128),  f{v,  ka)  = 
f{—i/,ka)).  The  residues  associated  with  the  poles  at  t/  =  «/„  are  given  by 


lim  (1/ -  i/„)/(i/,  ia)  = 

fco>l 


-1, 


(62) 


by  (153).  Thus  we  write 


f{i^,ka) 


7(1/,  ka)  ^  /  1 

ka  ;“j  V 1/  -  J/, 


1 


l'  +  I'n 
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(63) 


_  _  2  y' 

ka 

Following  similar  arguments  to  those  which  lead  to  (138)  we  find  that  use  of  this 
expression  gives  the  following  approximate  representation  of  ^*^(5,0): 

r{s,o)  v>"(s,o)  +  4f;i/„e“'’'i^i 

V>^(5,0)  +  -f;e-"l^l.  (64) 

“  n=l 

We  expect  this  model  to  be  reasonably  accurate  for  k\s\  ~  1.  Since  Imi^Vn)  >  0 
and  Im(un)  oc  (kaY^^  the  curvature  correction  term  is  strongly  attenuated  as  A:|s|  is 
increased.  Indeed,  upon  inspection  of  (58)  and  (64)  it  is  apparent  that  this  term  rep¬ 
resents  the  so-called  ‘missing’  creeping  waves  referred  to  earlier  in  this  section.  They 
are  numerically  significant  since  they  traverse  less  than  half  the  cylinder  circumfer¬ 
ence.  Although  we  expect  the  representation  (64)  to  be  satisfactory  for  ibls|  ~  1 
it  is  clear  that  it  breaks  down  for  k\s\  •C  1.  Indeed,  as  |s|  (or  equivalently  |0|  ) 
approaches  zero  the  imaginary  part  of  the  correction  term  blows  up,  a  feature  not 
present  in  plots  of  the  exact  form  of  ip‘{s,Q). 

It  is  worth  making  a  comment  here  on  the  analysis  of  4>’^{s,0)  given  in  Appendix 
B.  Although  the  split  into  two  components  and  (through  use  of  (130)) 
is  mathematically  convenient  it  is  not  unique  and  it  does  not  lead  to  separation  of 
physically  distinct  contributions  as  it  does  in  the  analysis  of  the  scattering  of  a  plane 
wave  by  a  circular  cylinder  [11].  In  particular  creeping  waves  appear  in  the  analysis 
of  both  components.  It  is  quite  possible  that  a  different  split  would  be  justifiable  on 
physical  grounds. 

So  far  in  this  section  two  simple  representations  ((57)  and  (64))  of  the  j/>-distribution 
of  the  circular  cylinder  have  been  obtained.  They  have  been  shown  to  be  valid  for 
different  values  of  A:|s|.  Ideally  we  would  like  a  single  expression  accurate  through¬ 
out  the  interval  0  <  A:|s|  <  wka.  Maystre  [1]  has  given  an  empirical  formula  which 
is  able  to  represent  ^’“’(•8,0)  to  a  very  high  degree  of  precision.  It  appears  that 
this  was  obtained  by  performing  a  seven  parameter  study  using  limited  information 
on  the  expected  anal}i,ic  form.  As  a  result  Maystre’s  formula  is  complicated  and 
may  not  be  readily  applicable  to  other  smooth  scatterers.  The  present  author  has 
attempted  a  rigorous  analysis  of  the  exact  form  of  V^'^(s,  0)  in  a  bid  to  discover  a 
simple  and  accurate  representation.  The  fact  that  this  approach  failed  to  produce  a 
single  usable  form  is  directly  traceable  to  the  absence  of  an  analytic  expression  for 
the  Fourier  transform  of  HY^'{z)I HY\z)  (with  respect  to  v)  accurate  for  all  real  v 
and  z.  Nevertheless,  as  we  shall  now  demonstrate,  it  is  possible  to  produce  a  very 
good  approximation  to  i/’‘’(s,0)  for  0  <  k\s\  <  nka  using  a  simple  single  parameter 
model.  This  is  achieved  by  making  use  of  much  of  the  information  gleaned  from  the 
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rigorous  analysis.  Specifically,  we  postulate 

(65) 

with 

VQ  =  ka^  ao(ita/2)’/V"/^  (66) 

The  correction  term  has  been  chosen  to  ensure  accuracy  for  A:|s|  =  fca|0|  ■C  1  and 
to  represent  the  combined  effect  of  the  creeping  waves  traversing  less  than  half  the 
cylinder  circumference.  When  the  parameter  ao  takes  the  value  0.55,  (65)  is  an 
accurate  representation  for  0  <  A:|3j  <  vka.  This  is  demonstrated  in  Fig. 8  for 
ka  =  5  and  A:a  =  10.  Fine  tuning  of  qq  may  lead  to  even  better  agreement  with  the 
exact  form.  Computation  of  the  cylinder  current  using  (65)  (which  we  shall  call  the 
improved  DBO  (IDBO)  approximation)  is  described  in  the  next  section. 
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6  Application  of  IDBO  Theory  to  Scattering  by 
a  Circular  Cylinder 

The  scattering  of  a  TM-polarised  plane  wave  by  a  perfectly  conducting  circular 
cylinder  has  been  examined  rigorously  in  Section  2  and  using  ADBO  theory  in 
Section  4.  The  currents  obtained  by  applying  the  ADBO  technique  were  found  to 
be  inaccurate  in  the  shadow  region.  Here  the  problem  is  revisited  using  the  imp  ''ved 
form  of  the  ^-distribution  for  the  circular  cylinder  obtained  in  Section  5.  From  (38) 
and  (65)  we  find 


jiDBoiO]  =  -i^cos^e-*'--*®-  ^/,(^,A-a)  -  ^A/,(^,A-a),  (67) 

with 

MA(d,ka)  = 

The  integral  can  be  analysed  using  asymptotic  techniques.  (See  .A-ppendix  C.) 
It  is  found  that 


Mr{e,ka) 


,co$  6^0 
ir” 


—  (ka)'^  sin^  0) 


K(e.ka)  +  Ar/{0.ka). 


(69) 


where 


and 


Miie.ka) 


{i/^  —  {ka^  sin^  0) 


(70) 


J\/;P(^,  ka)  L^(|^|-7r/2)2;r 


2_ 

ka 


1/3 


Ai(Qo£’"/")(e 


71) 


Since  j/q  ka  for  ka  \  and  using  also  the  result  (48)  for  jADBo{0)^  it  follows  that 


/n\  *a>l.co»e?iO 

JIDBo{0)  ~ 


-M  cos  (i  +  sn^) 

-fn(e,ka)-fMUO,kii) 

-•i(n{e.ka)  +  i«(eM)) 
-'■^(M‘,(e,ka)kM‘/(e,ka)) 


for  0  <  |0|  ~  T 12 


for  IT  12  ~  10|  <  ir 


(72) 


Several  comments  are  appropriate: 

1.  Comparing  (72)  w-ith  the  asymptotic  form  of  the  exact  current  jex{0)  (49)  we 
see  that  IDBO  theory  yields  the  correct  curvature  correction  to  the  PO  current 
in  the  illuminated  region.  Furthermore  the  erroneous  curvature  correction  term 
present  in  jADBoi^)  for  7r/2  ~  1^1  <  tt  is  precisely  cancelled  by  a  contribution 
of  similar  form  due  to  Mri0,ka). 
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2.  The  /*P  and  M’^  terms  together  generate  an  approximate  representation  of 

the  most  significant  creeping  wave  contributions  to  the  current  on  the  shad¬ 
owed  side  of  the  cyhnder.  In  Section  5  it  was  observed  that  creeping  waves 
traversing  more  than  half  the  cyhnder  circumference  were  severely  attenuated 
and  gave  ar  insignificant  numerical  contribution  to  the  DBO  T/>-distribution. 
Their  contribution  was  therefore  neglected.  Thus  no  creeping  wave  terms  are 
present  in  the  IDBO  prediction  of  the  current  in  the  lit  region.  In  Appendix 
D  it  is  confirmed  that  the  correct  creeping  wave  behaviour  for  |0|  ~  7r/2  is 
generated  by  retaining  the  term  0)  (58)  in  the  approximate  represen¬ 
tation  of  Likewise  all  higher  order  creeping  wave  contributions  to  the 

current  in  the  shadow  region  are  produced. 

3.  The  /‘  and  Ml  terms  are  present  because  the  approximation  (65)  provides 
an  inadequate  model  of  v'^(s,0)  for  s  ~  ±7ra.  This  is  not  a  serious  concern 
because  even  for  moderate  ka  they  are  numerically  very  small. 

It  has  already  been  confirmed  in  Section  4  that  the  substitution  of  for  in 
(41)  makes  little  numerical  difference.  In  a  similar  way,  and  in  particular  because 
Im(i’o)  >  0.  we  expect  that  the  replacement  of  Mr  by  M^  in  (67)  will  introduce 
little  error.  Thus  we  write 

JIDBO(O)'^^'  -p^(0.ka)-jM^{0,ka)  (73) 

and  we  may  employ  FFT  techniques  for  an  efficient  compulation  of  /oc  and  Moo-  In 
a  combined  evaluation  of  the  second  and  third  terms  on  the  right  side  of  (73)  use 
can  be  made  of  tiie  analytic  form  for  the  Fourier  transform  of  our  approximation 
(65)  to  vHs.O): 


y>(-)  =  r  v'[a0.Q)r'-'^d0 

*/  — OC 

ka>\  -)(w\A-a)  ^  k  i/Q 

*0^2  2^* 
a  Za  (i/^  — 

Note  that,  since  Im{i'o)  >  0,  the  second  term  is  not  singular  for  real  u'.  (Recalling 
that  (65)  was  derived  in  an  attempt  to  approximate  the  Fourier  transform  of  /(i',  ka) 
with  respect  to  f  we  expect  a  better  approximation  to  V’'(u-’)  to  be  given  by 


^’“^(0,’)  sr  —kf(ra,ka)  =  —k 


Hl'hka)' 


(75) 


However,  this  result  is  likely  to  be  less  readily  generalisable  for  application  to  an 
arbitrary  smooth  scatterer.) 

Plots  of  the  magnitude  and  phase  of  the  surface  current,  as  given  by  (67),  are 
shown  in  Fig. 9  for  ka  =  5  and  ka  =  10.  The  integrals  Jr  and  Mr  have  been  evaluated 


19 


numerically.  Also  shown  are  the  results  obtained  using  (73)  and  the  analytic  form 
(72).  As  expected  invoking  the  extra  approximations  h  «  loo  and  ~  Moo 
introduces  little  error.  It  is  also  clear  that  (72)  is  indeed  a  valid  representation  of 
the  cyhnder  current  obtained  using  the  IDBO  technique  away  from  the  region  of  the 
shadow  boundary. 

In  Fig. 10  the  cyhnder  currents  obtained  by  the  IDBO  method  are  compared  with 
the  corresponding  exact  and  PO  solutions.  In  general  the  IDBO  currents  are  much 
more  accurate  than  those  of  PO  at  all  points  on  the  cylinder  surface.  There  is, 
however,  some  residual  erroneous  oscillatory  variation  with  angle  in  the  shadowed 
region. 
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7  Prescription  of  the  IDBO  Technique  for  a  Gen¬ 
eral  Two-Dimensional  Smooth  Scatterer 


In  Section  6  the  IDBO  method  was  successfully  applied  to  the  problem  of  plane  wave 
scattering  by  a  perfectly  conducting  circular  cylinder.  This  should  not  have  surprised 
us  since  the  technique  was  based  upon  a  derivation  of  an  accurate  approximation 
to  the  known  V’'<Iistribution  for  that  scatterer.  The  technique  can  only  be  used  to 
advantage  when  its  validity  has  been  demonstrated  for  bodies  for  which  the  exact 
t/>-distribution  is  not  known.  Only  then  will  the  procedure  yield  a  true  predictive 
capability. 

The  simple  approximation  (65)  to  V’‘^(-s,0)  was  developed  with  the  aim  of  be¬ 
ing  readily  generalisable  for  application  to  an  arbitrary  2D  smooth  scatterer.  We 
therefore  postulate  the  following  form  for  the  t/>-distribution  of  such  a  body: 

i'lDBois^.S,)  =  V'^(Si,S2)  -  (76) 

4a(si,S2) 

where  oq  =  0.55  and  Si  and  sa  denote  the  arc  lengths  of  points  Pj  and  P2  on  the 
scatterer  surface  from  some  reference  point  also  on  that  surface.  This  expression 
is  sufficiently  general  to  allow  dependence  on  both  variables  sj  and  $2  rather  than 
just  the  single  quantity  si  —  $2-  We  have  yet  to  prescribe  the  functions  a(s),S2) 
and  ^(51,52).  Any  proposed  l:>rm  fo:  a(si,S2)  should  reduce  to  the  constant  value 
a  for  the  Ccise  of  a  circular  cylinder  of  radius  a  while  |if(si,S2)l  should  reduce  to 
{kl2a^YI^\s,-S2\. 

To  obtain  suitable  expressions  for  a(sj,  S2)  and  ^(s] ,  S2)  we  appeal  to  Fock  theory. 
In  the  application  of  Fock  theory  to  a  general  convex  body  [11]  it  is  seen  that  the 
amplitude  of  the  current  in  the  shadow  region  depends  only  on  the  radii  of  curvature 
of  the  surface  at  the  point  of  excitation  [P\  say)  of  the  creeping  wave  and  at  the 
point  of  observation  {P2)  and  also  on  the  quantity 


r>2 

i 


1/3 


ds'. 


(77) 


2a^s’)^ 

where  a(s)  is  the  radius  of  curvature  as  a  function  of  arc  length.  Thus,  in  order  also 
to  preserve  symmetry  under  interchange  of  Sj  and  $2,  we  propose  that 

a(-si,52)  =  ■^0(51)0(32),  (78) 

and 


({si,S2)  =  j^ 


k 


1/3 


ds'. 


(79) 


In  accordance  with  (37)  the  IDBO  prediction  of  the  surface  current  induced  by  an 
incident  field  E'  on  a  closed  body  of  perimeter  I  would  be  as  follows: 


d  E'  r-’+V^ 

jiDBois)  =  -^{s)  P  J  ipjDBo{s,s')E'(s')ds' . 


(80) 
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Thr  form  of  the  t/)-distribution  given  by  (76), (78)  and  (79)  has  yet  to  be  tested 
for  any  scatterer  other  than  the  circular  cylinder.  Suitable  test  targets  would  be 
the  elliptic  and  parabolic  cylinders  for  which  control  solutions  are  readily  available. 
(See,  for  example,  Ref.  12). 

For  a  general  smooth  body  a(si,S2)  and  ^(^ijSj)  will  be  complicated  functions 
of  both  5i  and  $2.  Thus  the  contribution  of  the  curvature  correction  term  of  xI^idbo 
to  the  current  jwBO  w’ill  not  be  amenable  to  rapid  approximate  calculation  using 
FFT  algorithms. 

The  prescribed  forms  (78)  and  (79)  ensure  that  T/’/r»Bo(-Si, sj)  is  symmetric  and 
so  the  fields  calculated  using  the  IDBO  scheme  are  guaranteed  to  comply  with  the 
reciprocity  principle.  (See  [1]  and  Section  3.) 
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8  Conclusions  and  Suggestions  for  Further  Work 

In  the  DBO  formulation  of  the  electromagnetic  scattering  problem  the  current  is 
expressed  in  the  form  of  an  integral  over  the  scatterer  surface  with  the  integrand 
containing  a  so-called  ^’'distribution  which  is  characteristic  of  the  scatterer  and  de¬ 
pends  only  on  the  polarisation  of  the  incident  radiation.  The  ^’-distribution  gives  the 
current  induced  at  a  point  on  the  scatterer  surface  due  to  a  delta-function  excitation 
at  a  second  such  point.  In  this  report  two  approximate  techniques  based  on  the 
DBO  formulation  have  been  examined.  These  schemes  aim  to  provide  accurate  rep¬ 
resentations  of  the  V’-distribution.  A  simple  approximation  (the  ADBO  procedure) 
is  to  use  the  known  ^’-distribution  for  the  infinite  plane.  Although  used  success¬ 
fully  in  the  past  to  treat  scattering  by  rough  planar  surfaces  [3,4]  we  saw  in  Section 
4  that  it  gives  unsatisfactory  predictions  of  current  and  bistatic  RCS  for  circular 
cylinders  in  the  resonance  region.  In  particular  it  includes  only  half  of  the  known 
curvature  correction  to  the  PO  current  in  the  illuminated  region  and  over-estimates 
tha  magnitude  of  the  shadow  region  currents.  An  improved  approximation  was  de¬ 
rived  and  applied  to  the  circular  cylinder  problem.  This  IDBO  scheme  led  to  good 
agreement  with  the  known  current  everywhere  on  the  cylinder.  A  generalisation  of 
the  IDBO  technique  to  allow  treatment  of  scattering  by  arbitrary  2D  smooth  bodies 
was  proposed  in  Section  7. 

In  general  the  approximate  DBO  techniques  have  the  following  features: 

1.  Unlike  PO  their  application  leads  to  predictions  of  non-vanishing  surface  cur¬ 
rents  in  unlit  regions  and  smooth  transition  currents  at  shadow  boundaries. 

2.  Unlike  PO  they  yield  currents  which  do  not  depend  solely  on  the  local  geom¬ 
etry  of  the  scatterer;  i.e.  the  currents  are  non-local.  They  make  explicit  the 
phenomenon  of  short-range  coupling,  demonstrating  that  the  surface  current 
at  any  point  on  the  body  is  influenced  by  the  shape  of  the  scatterer  within  a 
distance  along  the  surface  of  about  a  wavelength. 

3.  A  single  expression  gives  the  surface  current  everywhere  on  the  body.  This  con¬ 
trasts  with  Fock  s  calculation  of  surface  currents  in  the  vicinity  of  the  shadow 
boundary  [11].  Interpolation  formulae  are  required  to  extend  Fock’s  results 
into  the  deep  illumination  region. 

4.  Provided  that  the  assumed  ^’-distribution  is  symmetric  with  respect  to  its 
arguments  the  calculated  fields  will  comply  with  the  reciprocity  principle. 

5.  They  are  more  efficient  than,  for  example,  integral  equation  techniques  because 
the  currents  are  given  explicitly. 

6.  They  are  readily  applicable  for  any  incident  field. 

Any  attempt  to  produce  accurate  current  and  field  predictions  using  approximate 
DBO  techniques  must  be  based  upon  a  reliable  representa'.  ion  of  the  appropriate 
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^-distribution.  It  seems  likely  that  an  accurate  model  of  the  ^’-distribution  for  a 
general  scatterer  many  wavelengths  in  size  can  be  generated  by  gleaning  information 
from  the  V’-distributions  of  a  small  number  of  standard  (canonical)  bodies.  This  is 
the  philosophy  of  the  successful  Geometrical  Theory  of  Diffraction  and  is  also  the 
reasoning  behind  the  in-depth  investigation  of  the  ^’-distribution  of  the  circular 
cylinder  in  this  paper.  An  IDBO  treatment  of  scattering  by  bodies  with  edges 
requires  an  examination  of  the  ^’-distribution  of  an  infinite  wedge.  This  work  is  well 
advanced  and  will  be  reported  in  a  future  paper.  An  alternative  approach  would 
be  the  approximate  numerical  computation  of  the  i/>-distribution  of  an  arbitrary 
scatterer  by  disregarding  the  (usually  insignificant)  influence  of  distant  parts  of  the 
body. 

In  this  paper  we  have  been  concerned  only  with  the  case  of  TM  incidence  on  a 
perfectly  conducting  infinite  cylindrical  body.  Provided  that  the  IDBO  technique  is 
properly  validated  for  this  case  it  seems  that  further  development  of  the  technique 
to  allow  tT-eiilment  of  TE  incidence,  full  3D  vector  scattering  and  possibly  scattering 
by  dielectric  and  magnetic  objects  may  be  both  feasible  and  beneficial. 
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A  Analysis  of  the  Integral  7;r(^,  ka) 


The  integral  /,(^,  jka)  (Eq.(42))  occurs  in  the  ADBO  treatment  of  scattering  by  a 
conducting  circular  cylinder.  Exact  analytic  evaluation  of  this  integral  appears  not 
to  be  possible.  Here  we  attempt  an  approximate  ev^aluation  for  A:a  ^  1. 

The  integrand  of  1^(0 ,  ka)  behaves  as  (1/^')^  for  C  1  and  so  it  is  not  integrable 
in  the  usual  sense.  However,  as  discussed  in  detail  by  Maystre  [3],  integrals  of  this 
type  arise  in  the  DBO  technique  as  derivatives  of  principal  value  integrals  and  so 
they  are  well-defined.  Thus  we  write 


ln{0,  ka)  =  lim 


r-o  Trfcae 


^—ikacosff  _j_ 


(/  Ic)  ■  (81) 


To  render  1^(6,  ka)  in  a  form  suitable  for  asymptotic  analytic  evaluation  (and  indeed 
numerical  evaluation)  w'e  use  Bessel’s  equation  to  give 

I,(e,ka)  =  i  ^  +  {kafHi,'\kae'))  (82) 

and  then  integrate  by  parts  to  find 

7,(0,  ka)  =  i  f  cos(0'  -I-  0)(1  -  ika  cos(0'  +  6))de' 

-  277|^^(7rita)e“""°*®.  (83) 

The  integrand  of  the  remaining  integral  possesses  a  logarithmic  singularity  at  0'  =  0 
and  so  is  integrable  in  the  conventional  sense.  Asymptotic  analysis  of  this  integral 
is  most  convenient  through  consideration  of  the  integral  /v„(0,/:a): 

K,{e,ka)  =  I'  (84) 

Straightforward  manipulation  shows  that  /»{0,  ka)  may  be  expressed  in  terms  of 
K-ir{0ika)  as  follows; 

J,{e,ka)  =  /'-(«,  to).  (85) 

When  ka  >  1  the  function  Ho^kaO')  is  not  slowly  varying  in  the  region  -tt  < 
6'  <  TT.  However,  asymptotic  evaluation  of  Kr{6^ka)  can  proceed  by  conventional 
techniques  using  the  following  integral  representation  [13]: 

Hi^\kae')  ^  —  r  (86) 

ITT  J — oo 

Substituting  into  (84)  and  interchanging  the  order  of  integration  gives 

h%{0,ka)  =  l<l-^\0,ka)  +  Ki~H0^ka).  (87) 
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where 


KMie,ka)  = 

ITT  J-oo  Jo 

(88) 

Kl-He,kc)  = 

1  foo  fO 

T-  /  / 

ITT  J— oo  J— IT 

(89) 

±0'cosht  —  cos{0'  -f  0). 

(90) 

A  proper  asymptotic  analysis  of  these  double  integrals  must  account  for  both  sta¬ 
tionary  phase  point  and  end-point  contributions.  (See,  for  example,  Ref.l4.)  These 
contributions  are  now  considered  in  turn. 

A.l  Stationary  Phase  Point  Contributions 

For  evaluation  of  the  stationary  phase  point  contributions  to  Ktt{6,  ka)  it  is  conve¬ 
nient  to  note  that  is  even  in  6  ;  i.e.  ka)  =  !{■„{— 6,  ka),  and  therefore  that 

the  analysis  can  be  confined  to  a  treatment  of  the  interval  0  <  0  <  tt.  Consider  first 
the  double  integral  ka).  For  a  stationary  phase  point  at  (ff\i)  =  (^o,<o)  the 

following  pair  of  simultaneous  equations  must  be  satisfied: 

cosh  to -f  sin(0o  4- 0)  =  0  (91) 

0osinhto  =  O.  (92) 

A  solution  is  possible  only  when  to  =  0  and  sin(^o  -(■  ^)  =  — 1-  It  follows  that  there 
is  no  stationary  phase  point  contribution  to  ka)  in  the  interval  0  <  &  <  7r/2. 

For  7r/2  <  ^  <  TT,  however,  there  is  a  stationary  phase  point  in  the  integration 
interval  0  <  ^'  <  rr  at  =  37r/2  —  0.  To  evaluate  its  contribution,  which  we  shall 
call  ka),  we  introduce  a  new  integration  variable 

0"  =  0' -0'^  =  0'  +  0-  37r/2  (93) 

so  that 

1  foo  eB—rl2 

Ki^^0,ka)=  I  I  '^U9"dt.  (94) 

ITT  J~<x  Je~3ir/2 

Expanding  the  exponent  in  a  Taylor  series  about  the  second  order  stationary  phase 
point  at  {0",t)  =  (0,0)  we  find 

9^(0",  t)  «  (37r/2  -  ^)(1  -i-  tV2)  +  ^'76  +  0''ty2.  (95) 

The  t-integral  can  be  performed  using  the  standard  result  [14] 
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Thus 


),,«  t'-’n  .fg" 


A-<+>-<'(.«,J:a) 


/2_  ;.e-5r/2  gisaff  /D  jun 

- /  ,  for  7r/2  <  0  <  tt. 

ITT  V  Kfl  Je-3>r/2  +  37r/2  -  0 


V*'  ■  ! 

For  Ara  1  the  leading  contribution  to  the  remaining  integral  is  of  order 
This  is  the  contribution  due  to  the  stationary  phase  point  at  6"  =  0: 

I  ^  ^  ^  I  ^ikaB'^le 

•'«-3»/2  j0"  +  37r/2-0  ~  ^37r/2  -  0 


/ZTTl2-eJ-°^ 

2ir  /  2 

,  [Pi  ^^(0). 

rx-rr  0  -  ft  \kaj 


y37r/2-0 

Here  we  have  used  an  integral  representation  of  the  Airy  function  Ai{x)  [14].  The 
constant  Ai(0)  is  approximately  equal  to  0.355.  The  effects  of  the  finite  integration 
limits  will  be  considered  more  precisely  (i.e.  without  the  stationary  phase  approxi¬ 
mation  to  the  integrand)  later  in  this  Appendix.  We  have  found  that 


for  0  <  0  <  7r/2 


3)T/2-e 


A  similar  determination  of  the  stationary  phase  point  contribution  to  the  integral 
ka)  is  possible.  Again  a  stationary  phase  point  exists  within  the  integration 
interval  only  when  7r/2  <  0  <  tt.  It  is  located  at  {0\t)  =  (7r/2  —  0,0).  The  change 
of  variable  6"  —  6'  6  —  n /2  gives 


1  foo  r$—'n/2 

Ki-\e,ka)  =  —f  / 

ITT  J^ooJo~3n/2 


oo  3ir/2 


(100) 


Using  (90),  expanding  about  the  stationary  point  at  {6",t)  =  (0,0)  yields 

g.{e'\i]  %  (0  -  7r/2)(l  -I-  fV2)  -  0'76  -  e"t^l2.  (101) 

Carrying  out  the  f-integration  as  before,  we  find 


Kl-'-'ie.ta) 


k.>l  £-'*■»  /«  M" 


llni^  /  . .  for  nl2<6<  n. 

ka  ye-37r/2  _  7r/2  -  0" 


(102) 

This  time,  however,  the  function  multiplying  the  exponential  in  the  integrand  is  not 
well-behaved  near  the  stationary  phase  point  0"  —  0  when  6  k  n  12.  Thus 


r«-V2  dO"  ka>i  2n 


-3^/2  J0-T^f2-  0" 


27!-  /  2 

-  (  —  )  ^J(O) 

le-  7r/2 


(103) 
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provided  ^  96  7r/2.  Hence 


fca>l,«94ir/2 


for  0  <  6  <  7r/2 


[  (li)  fe)  '  for»/2<«<. 

(104) 

and  so  from  (87), (99)  and  (104)  we  deduce  that 


K:’{0,  ka)  t,(|s|  _  ^/2) 


g.fca(34r/2-|fi|)  g.fca(ie|-4r/2) 

JiTT/i  -  1^1  J\0\  -  ,1-1 


(105) 


Substituting  into  (85),  we  find  the  following  contribution  to  I^{6,ka)  arising  from 
stationary  phase  considerations: 


inSM)  -  ,/2)  (1^)''"  (^)’'“/li(0)e-=‘’l- 


g«A:a(3)f/2-|e|)  g»7co(|e|-»r/2)  \ 

(3i7/2  -  (191  -  »/2)3/2  j  ■ 


(106) 


A. 2  End-Point  Contributions 


Asymptotic  series  for  the  end-point  contributions  to  the  6'  integrals  of  (88)  and  (89) 
can  be  developed  using  successive  integrations  by  parts: 


r  ^  f  ,  0  fiy  1 

Jo  ikag'^{7r,t)\  ika{g'^{n,i)y  \ka)  j 

,.fca,,(0,t)  r  I  o(^  M'] 

ikag\.{0,t)  \  ika{g\.{0,t)y  \ka)  J 

r  (1  +  +  0  f  i-Vl 

J-TT  *^0^1(0, t)  \  ika{g'_{0^t))'^  \kaj  J 


(107) 


I,  +  8-(-».0  ,  0  I'-V'l  .(1081 

ikag'_{—'n',t)\  ika{g'_{—'ir,t)y  \/:a/  J 


Here  g'^  and  g'^  respectively  denote  the  first  and  second  derivatives  of  g±  with  respect 
to  0'.  Hence  the  end-point  contributions  to  K-n{0.,ka)  (denoted  by  K‘{d,ka))  are 
found  to  be 


K{6M)  = 


6-**°*=°*^  f  p  /  1  1  \ 

Tcka  \*/— OO  \cosht4-sin0  cosh  <  — sin  0/ 
COS0  f  ^  ^  ^ 

ika  J-oo  \  (cosh  t  +  sin  6)^  (cosh  t  —  sin  0)^  / 

r  /"  (  1  ^  1  N 

irka  U-oo  \ cosh <  + sin 0  cosht  — sin0/ 

cosd  /■oo  /  1  1  \ 

ika  j-oo  \ (cosh  i  +  sin  0)^  (cosh  <  — sin  0)^/ 

+  “(s)  }■ 


sin  6)' 


pikaco»e  (  ,oo 


dt+0 


irka  cosht. 


-irka  coiht. 


(109) 


The  first  term  is  due  to  the  discontinuity  of  the  derivative  of  the  integrand  of 
Kr{0,ka)  at  6'  =  0.  The  first  ^-integration  can  be  performed  by  application  of 
contour  methods  in  the  complex  plane  or  by  repeated  use  of  the  standard  integral 
[15] 

dt  a  .  . 


cosh  t  +  cos  a  sin  a 


(0  <  Q  <  tt). 


(110) 


We  find 


/*=  /  1  1  \  27r 

-oo  \ cosh  f -1- sin  0  cosh  f  —  sin  0/  lcos01  ^  ^ 

Differentiating  this  result  twice  with  respect  to  sin  0  gives 

r(  1  I  1  ^  (3-2cos^0) 

V-oo  \(cosh  f -I- sin  0)^  (cosh  f  —  sin  0)^y  |  cos  01  cos'*  0 

The  second  term  of  (109)  is  due  to  the  truncation  of  the  integrand  at  the  end-points 
0'  =  ±7r.  Provi.  ed  |sin0|  76  1  the  ^-integrals  can  be  evaluated  for  ka  I  by  the 
method  of  stationary  phase.  Expanding  the  integrands  around  the  stationary  phase 
point  at  t  =  0  in  the  usual  way  [14]  we  find 

/I  1  \  .Vl-n  *  I.  /cO  >  1  .COS  ffStO 


r  ( _ u 

J  — oo  \(cosh  t  +  sii 


(111) 


(112) 


r  ( _ L_ 

J  —oo  Vcosh  f  -f  si 


- V -  e 

sin  0  cosh  i  —  sin  0/ 


•TrJba  cosh  t  Jj 


—  _ I 

ka  cos 


1  f  1  _  t(8-3cos^0) 
cos^  0  I  Stt/co  cos^  0 


/j_V 

\/:a/ 


(113) 


/:( 


(cosh  < -f- sin  0)^  (cosh  f  —  sin  0)' 


^irkacosht^^  ka^l^sB^O 


(As  before,  (114)  can  also  be  obtained  by  differentiating  (113)  twice  with  respect 
to  sin0.)  Substituting  (1 1 1  ),(1 12),(1 13)  and  (114)  into  (109)  gives  the  end-point 
contribution  to  Kr{0,ka): 
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/cal  COS  I  2tkacos^6 


ih] 


/  2  ^ika  coi  8  ^inka^iir/ 4 

V/ca/  IT  COS^  0 

j  i(8  — 3cos^0)  z(4  — 3cos^0) 


Sirka  cos^  0 


ka  cos^  0 


We  now  substitute  this  result  into  (85)  to  find  the  end-point  contributions  to  lr{0,  ka) 
(denoted  /*(^, /ca)).  After  performing  the  necessary  differentiations  and  much  te¬ 
dious  algebra  we  find 


l‘{0,ka)  2|cosf?|e 


— f^o  COS  B 


2ka  cos^  0 


+/;(o,*-a)  +  o(i)\ 


(116) 


TV  i  tii 

(0  \  3/2  pZi-n / 4 piirka  pika  coi  6 

We  have  included  the  first  term  on  the  right  side  of  (85)  and  have  used  the  asymp¬ 
totic  form  of  the  Hankel  function  to  arrive  at  this  result.  This  is  appropriate  since 
that  term  originated  from  end-point  contributions  when  integrating  by  parts.  It  is 
of  interest  to  note  that  the  term  7^(0, /ca),  which  is  due  to  the  truncation  of  the 
integration  interval  at  0'  =  ±7r,  can  also  be  derived  directly  by  considering  loo  —h- 

loo{0,ka)-h{0,ka)=^  (^J~\  /")  (118) 

Using  the  asymptotic  form  of  the  Hankel  function  we  find 

I^{e,ka)-h{e,ka)  fee’""'*/ 


Thus 


as  required. 


\lTrka  Jr  {0'fl^" 

co._£5i0  /  2  ^-,r/4^,rka^,kacose 

V/ca/  TT^COS^^ 

l,{0,ka)'''^'  loo{0,ka)  +  n{0,  ka) 


[if- 


(119) 


(120) 


Finally,  then,  we  add  the  stationary  phase  and  end-point  contributions  to  deter¬ 
mine  the  asymptotic  approximation  to  l„{0,ka): 


7,(0, /ca)  r/{O,ka)^It{0,ka). 
where  P/  and  7'  are  given  by  (106)  and  (116)  respectively. 


(121) 
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B  Asymptotic  (ka  >  1)  Analysis  of  'ip^(s^O) 


Our  starting  point  is  the  exact  expression  (32); 


^=(s,0) 


-k 

2x0 


OO 


(122) 


with 

6  =  s/a 


(123) 


and 


ni.kc) 


(124) 


An  asymptotic  investigation  is  possible  by  means  of  a  Watson  transformation  [II]. 
Using  this  technique  the  infinite  sum  is  converted  to  a  contour  integral  in  the  complex 
plane  (see  Fig.  11): 


E 


/jS-OC 


L 


C+  +  C- 


sin  VTT 


/(i2,  ka)di/. 


(125) 


This  step  is  e£isily  checked  by  evaluating  the  residues  of  the  (simple)  poles  of  sin  i/r 
which  lie  on  the  real  axis.  The  assumption  has  been  made  that  f{u,ka)  has  no 
poles  lying  within  the  contour  C+  +  C_.  By  choosing  the  positive  constant  a  to  be 
arbitrarily  small  this  is  indeed  the  case  —  Hl^^ka)  =  0  has  no  solution  for  real  v. 
Since  [14] 


H»\ku)  = 

(126) 

//!"'( to)  = 

(127) 

we  deduce  that 

fy.ka] 

=  /(-I/,  ka). 

(128) 

Hence 

cos.(«+o) 

sm  i/x 

(129) 

Following  the  analogous  treatment  of  the  field  scattered  by  a  perfectly  conducting 
cylinder  due  to  plane  wave  incidence  [11]  we  choose  to  write 


cos  t'(^  +  x)  =  — le  "'^sin  i/x  +  e'*'’"  cos  (130) 


and  examine  the  contributions  of  these  terms  separately.  It  remains  to  be  seen 
whether  similar  physical  interpretations  can  be  given  to  these  contributions. 

The  integrand  of  (^29)  due  to  the  first  term  on  the  right  side  of  (130)  has  no 
poles  lying  on  the  real  axis.  Thus  the  integration  contour  may  be  taken  to  be  that 
axis: 

v^<''(^,0)  ~  r  f{y.ka)du.  (131) 

^  TT  O  J  oc 
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As  far  as  the  author  is  aware  the  Fourier  transform  of  f{u,ka)  is  not  known.  To 
proceed  further  we  thus  revert  to  a  uniform  asymptotic  expansion  of  f{i/,ka).  Ap¬ 
plying  known  uniform  asymptotic  expansions  of  both  Hl^^'{ka)  and  Hl^^{kc)  [16],  it 
follows  that 


_  Hl}^'{ka)  Aro>i,li-l~i:o  '){v,ka)  ^  ka 

~  Hl'\ka)  ~  ^  -  {kayy 


(132) 


where  7(1/,  A:a)  is  defined  in  (19).  The  error  is  of  order  0(|i/|)  ^  or  0{ka)  Inserting 
the  leading  order  term  into  (131)  we  obtain  a  contribution  to  0)  of 


1 

27ra^ 


e~'‘"*/“7(i/,  ka)di/ 


-ikH[^\ks) 

~Y  s 


(133) 


where  the  integral  has  been  evaluated  using  (23)  and  (24).  All  dependence  on  ka 
has  disappeared  and  w'e  have  recovered  t/’'’(s,0).  The  next-to-leading  order  term  of 
(132)  contributes 

-P 

17  i-oc  {v^-(kar^r- 

To  evaluate  this  integral  we  need  to  know  how  to  treat  the  poles  at  1/  =  ±ka. 
This  problem  is  resolved  by  considering  their  origin.  They  have  arisen  because  we 
have  used  the  approximation  (132).  However  a  closer  examination  of  the  zeroes 
Un{n  =  1,2, . . .)  of  Hl^'>{ka)  [17]  reveals  that  (for  /?e(i^)  >  0)  they  lie  approximately 
along  a  line  oriented  at  60*’  to  the  positive  real  axis: 

i/n  =  ’-a  +  One’"/"  j  +0(ka)-^/\  (135) 


The  real  constants  Qn  are  related  to  zeroes  of  the  Airy  function: 


A2(-q„)  ==  0. 


(136) 


As  the  Q„  are  positive,  /m(//„)  >  9.  Thus  the  pole  cf  (1 34)  at  1/  =  ka  should  be  taken 
to  lie  just  above  the  real  axis.  Also,  in  view  of  (128),  the  pole  at  1/  =  —ka  is  taken 
to  He  just  below  tiie  real  axis.  Evaluation  of  (134)  then  follows  by  straightforward 
application  of  Cauchy’s  integral  formula: 


/ 


:di/  = 


—ik 


47r  J-oc  (i^^  —  (kaY)  4a 
’’om  (131),  (132),  (133)  and  (137),  then,  n  follows  that 

^’^<*'(5,0)  =  C’’(s.0)  -  +  k\0{ka)-\ 

4a 


(137) 


(138) 


The  second  term  represents  a  first  order  curvature  correction  to  v’^(^,0). 
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We  now  consider  the  contribution  of  the  second  term  of  (130): 

(139) 

The  integrand  has  simple  poles  at  the  zeroes  of  sini/7r  and  at  i/  =  where 

=  0.  Referring  to  Fig. 12,  we  see  that 


i 


cos  t/6 
(1  -  e2“'") 


/(i/,  ka)du  =  0, 


(140) 


where  the  contour  C  is  given  by  C  =  C+  +  Ci  +  C2  +  C3.  Consider  the  contributions 
of  the  arcs  at  infinity  (Ci  and  C3).  From  (132), 


f{y,ka)  ~  -j—.  (141) 

Since  Im{v)  >  0  and  |0|  <  tt,  this  behaviour  is  swamped  by  a  decaying  exponential 
due  to  the  rest  of  the  integrand.  It  is  therefore  clear  that  the  integrals  along  C\  and 
C3  vanish.  Thus 


^o(2)(^^0)=—  / 

Tra  Jc 


e^“"'  cos  v9 


c,  (1 


021m 


f{u,  ka)du. 


(142) 


Applying  Cauchy’s  integral  formula  once  again  we  find 

^^(2)(,,0) 


(143) 


where  the  sum  is  over  the  residues  of  the  poles  at  A  Wronskian  relation  for 
Hankel  functions  [17]  enables  us  to  write 


Hlinka) 


4i 


nkaHlV{ka)lHl^\ka)l=,„ 

James  [14]  gives  the  following  useful  asymptotic  forms  for  Hankel  functions: 

2 


(144) 


Hl^\ka) 


Mai 


(145) 

(146) 


From  (145)  we  deduce 


dv 


Hl^\ka)  2  -  ka)  (^Ay^%2.r/3| 


(147) 


Since,  for  large  fca,  w  ka  and 


/  9  \  1/3 

(148) 

(see  (135)),  we  find 

Hil\ka)  ‘“ft'  2  ^^y^%‘^Mi(Q„e-’^) 

(149) 

(150) 

Thus,  from  (144), 

//(i)'(/:a)  ^a>i  i 

f;:Hi^\ka)U,„  ''  27re2W3yij(o„e— /3)A2'(-q„)’ 

(151) 

Invoking  the  Wronskian  relation  [18] 

x)-^  Ai{xt'^l^)Ai\ 

27r 

(152) 

and  (136)  this  becomes 

H't}J-(ka)  j 

(153) 

Substituting  into  (143), 

“S'  2=*-f; 

a  (l-e^*""') 

(154) 

Since  >  0  this  may  be  rewritten  as 

^(e*‘'n(2mx+e)  (255^ 

^  n  =  l  m  =  l 

This  form  clearly  suggests  a  creeping  wave  interpretation.  Recalling  that  ss  ka 
for  large  ka  we  see  that  when  0  =  0  the  phase  factors  correspond  to  waves  which 
have  completely  encircled  the  cylinder  m  times  in  both  a  clockwise  and  anticlockwise 
sense. 

Finally,  then,  from  (138)  and  (155)  it  follows  that 
V-'(s,0)  =  + 

T^''’(s,0)- — 

4a 

°  n=l  m  =  l 
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C  Analysis  of  the  Integral  MT:{6^ka) 

The  integral  M^iO,  ka)  (Eq.{68))  occurs  in  the  IDBO  treatment  of  scattering  by  a 
conducting  circular  cylinder.  Exact  analytic  evaluation  appears  not  to  be  possible. 
Here  we  give  an  approximate  evaluation  valid  when  fca  1.  We  write 


Mr{9,ka)  =  Mi-^\9,ka)  +  Mlr\9M). 

(157) 

Jo 

(158) 

Ml-\9,ka)  =  /° 

J  —It 

(159) 

h±{9')  =  ±i^£)0'-icos(0'-he). 

(160) 

Since  M-n{d,  ka)  =  A/,(— ka)  we  confine  the  analysis  to  the  Ccise  0  <  0  <  tt.  Saddle 
point  and  end-point  contributions  are  now  considered  in  turn. 

C.l  Saddle  Point  Contributions 

Consider  first  the  integral  M)^'^\6,ka).  For  a  saddle  point  at  9'  —  6q  (possibly 
complex)  we  have 

^+(^o)  =  ^  (^)  +  ism{0'Q  6)  =  0.  (161) 

Since  i/Q  «  ka  for  /:a  ;§>  1  (see  (66))  solutions  to  this  equation  occur  in  closely  spaced 
pairs  in  the  complex  6'  plane.  It  is  not  permissible  to  treat  them  separately.  Instead, 
following  arguments  similar  to  Jones  [19],  we  expand  about  their  mid-point  9'  =  9[, 
defined  by 

/t"(0;)  =  tcos(<?;  +  0)  =  0.  (162) 

For  0  <  ^  <  7r/2  this  equation  has  no  solution  in  the  integration  interval  0  <  9'  <  z 
close  to  a  saddle  point  pair.  For  z  12  <  9  <  z^  however,  there  is  an  appropriate 
solution  at  9\  =  Zz {2  —  9.  The  associated  saddle  points  P±  are  at 

9'  (163) 

and  the  steepest  descents  path  passes  through  P±  with  slopes  tan"^(7r/6  ^  7r/4). 
Thus,  for  7r/2  <  9  <z, 

J  c 

where  t  =  9'  —  9[  and  the  contour  C  =  •+  -f  Ci  +  -f  .  (See  Fig. 13.1  The 

contributions  from  the  arcs  at  infinity  (C^  and  C^)  vanish.  The  contributions  from 
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Ct  and  Ct  represent  the  effects  of  the  truncation  of  the  integrand  of  ka) 

at  6'  =  ±n  and  will  be  considered  more  precisely  (i.e.  without  the  saddle  point 
approximation  to  the  integrand)  later  in  this  Appendix.  Deforming  the  steepest 
descents  contour  Ci  and  using  the  substitution 


z  = 


(165) 


yields  the  saddle  point  contribution  to  ka)\ 


jba)  U{e  ~  7r/2)(-i)  ^Ay^%'-o(3^/2- 


B) 


f  +/ 

Joce~'^/^  Jo 


roc  ' 


(166) 


The  remaining  integrals  furnish  a  familiar  representation  of  the  Airy  function  [14]. 
We  find  (for  0  <  ^  <  tt) 


Mi-^^‘^{e,ka)  U{0-7r/2)  27r 


An  analogous  treatment  of  the  saddle  point  contributions  to  Ml~{6,ka)  is  pos¬ 
sible.  The  result  (again  for  0  <  0  <  tt)  is  given  below: 

Mi-^’^{e,ka)^"^^  U{e-nl2)  27r  (168) 


C.2  End-Point  Contributions 


Asymptotic  series  for  the  end-point  contributions  to  the  integral  Mr{6,ka)  can  be 
developed  in  the  usual  way  (see  Appendix  A).  We  find 


M‘{e,ka] 


*0>1 


2ii/o 


—  {kay  sin^  6) 


(e 


—ika cos  $ 


—  e 


ivo'”  „ikacos0 


)• 


(169) 


The  first  term  arises  from  the  discontinuity  in  the  derivative  of  the  integrand  at  ^'  =  0 
while  the  second  term  is  due  to  the  truncation  of  the  integrand  at  the  end-points 


0'  =  ±7r. 


Finally,  then,  we  add  the  saddle  point  and  end-point  contributions  to  determine 
the  asymptotic  approximation  to  M^{0,ka): 

A/,(0,Jfca)  Mi-^^’^{0,ka)  +  Ml-^’^{0,ka)  +  M‘{0.ka),  (170) 

where  and  Af'  are  given  by  (167),  (168)  and  (169)  respectively. 
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D  Contribution  of  to  the  Cylinder  Cur¬ 

rent 

Here  we  compute  the  DBO  contributions  to  the  cylinder  current  due  to  creeping 
waves  which  traverse  more  than  half  the  cylinder  circumference.  Such  waves  arc 
accounted  for  by  the  term  of  (57).  We  begin  with  the  form  (154): 


,c(2)/  ^  COS  i/„0 

.  S  (1-eW) 

(171) 

From  (38)  we  deduce  that  this  generates  a  contribution  to  the  current: 

/“!(«) ‘I?'  A 

(172) 

where 

N,{0,ka)  =  /’(e"-"*'  + 

J  —  IT 

(173) 

This  integral  bears  many  similarities  to  the  integral  A/t(0,  ^ol)  analysed  in  Appendix 

C.  An  analogous  asymptotic  treatment  is  possible.  We  find 

NA0,  ka)  A7(0,  ka)  +  K(0,  ka), 

(174) 

where  the  end-point  and  saddle  point  contributions  are  given  respectively  by 

ATt  / /)  7.  _\  4t'n  sin  _»fcacos9 

A,r(,Vi  ~  /  0  /;  •  0\2\^ 

(i/2  -  {ka  smOy) 

(175) 

and 

A7(0,Jtfl)  - 

,  27r  Ai(Q„e*^)(e“'"<l«l-’7/2)  -|-  e»'"(3»/2-l»l)) 

for  0  <  |01  ~  7r/2 

for  7r/2  ~  |0|  <  T 

(176) 

The  end-point  term  arises,  as  usual,  because  of  the  truncation  of  the  integration 
interval  at  6'  =  ±n  and  contributes 


j/„e* 


-  (^'fl  sin  0)2) 


(177) 


to  the  current.  As  expected  the  magnitude  of  this  contribution  is  negligibly  small. 
The  contributions  of  the  saddle  point  terms  (176)  to  the  cylinder  current  are 


(178) 
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where  we  have  used  the  relation 

Ai{ane"l^)Ai\-ar.)  =  (179) 

(This  follows  from  (152)  using  Ai{—an)  —  0.)  Comparing  with  the  asymptotic  form 
(49)  of  the  exact  current  we  see  that  we  have  generated  the  correct  creeping  wave 
series  in  the  illuminated  region.  Furthermore,  by  using  the  expansion 

(1  _  2 

m=0 

(valid  since  /m(i/„)  >  0)  it  is  straightforward  to  confirm  that  the  precise  form  of 
all  higher  order  (i.e.  m  >  0)  creeping  wave  contributions  to  the  current  on  the 
shadowed  side  of  the  cylinder  has  been  reproduced  in  (178).  These  terms  represent 
waves  which  have  completely  encircled  the  cylinder  at  least  once. 
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Fig.l.  Electromagnetic  scattering  problem:  infinite  cylinder  illuminated  by  plane  wave. 


Fig. 2.  DBO  problem:  infinite  cylinder  illuminated  by  fictitious  source. 
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Fig.4(a).  ADBO  predictions  of  circular  cylinder  current  {ka  =  5). 

-  using  7,  (eq. (41)) 

.  using  7oo 

-  -  -  -  using  analytic  approximation  to  h  (eq.(48)) 
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Fig. 4(b).  ADBO  predictions  of  circular  cylinder  current  {ka  =  10) 

-  using /,  (eq. (41)) 

.  using  loo 

-  using  analj'tic  approximation  to  7^  (eq.(48)) 


X:  Blstatlc  Angle  Cdeg.D 
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Fig. 9(b).  IDBO  predictions  of  circular  cylinder  current  (ka  =  10) 

-  using  h  and  (eq.(67)) 

.  using  loc  and  (eq.(73)) 

-  -  -  -  using  analytic  approximation  to  h  and  (eq.(72)) 
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Fig. 11.  Integration  contour  for  Watson  transformation. 


Fig. 12.  Integration  contour  C  of  Eq.(140). 


t-plane 


contour  C  of  Eq.(164). 
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